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ABSTRACT 

We develop a new ID photochemical kinetics code to address stratospheric 
chemistry and stratospheric heating in hot Jupiters. Here we address optically 
active S-containing species and C0 2 at 1200 < T < 2000 K. HS (mercapto) 
and S 2 are highly reactive species that are generated photochemically and ther- 
mochemically from H 2 S with peak abundances between 1-10 mbar. S 2 absorbs 
UV between 240 and 340 nm and is optically thick for metallicities [S/H] > 
at T > 1200 K. HS is probably more important than S 2 , as it is generally more 
abundant than S 2 under hot Jupiter conditions and it absorbs at somewhat redder 
wavelengths. We use molecular theory to compute an HS absorption spectrum 
from sparse available data and find that HS should absorb strongly between 300 
and 460 nm, with absorption at the longer wavelengths being temperature sen- 
sitive. When the two absorbers are combined, radiative heating (per kg of gas) 
peaks at 100 /xbars, with a total stratospheric heating of ~8 x 10 4 W/m 2 for a 
jovian planet orbiting a solar-twin at 0.032 AU. Total heating is insensitive to 
metallicity. The C0 2 mixing ratio is a well-behaved quadratic function of metal- 
licity, ranging from 1.6 x 10~ 8 to 1.6 x 10~ 4 for -0.3 < [M/H] < 1.7. C0 2 is 
insensitive to insolation, vertical mixing, temperature (1200 < T < 2000), and 
gravity. The photochemical calculations confirm that C0 2 should prove a useful 
probe of planetary metallicity. 

Subject headings: planetary systems — stars: individual(HD 209458, HD 149026) 
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1. Introduction 

Stratospheric temperature inversions are ubiquitous in the Solar System, and it 
is beginning to look as if they are commonplace on hot Jupiters as well. Stratospheric 
temperature inversions form when substantial amounts of light are absorbed at low pressures 
(high altitudes) where radiative cooling is inefficient. Hubeny et al. (2003) pointed out that 
efficient absorption of visible light by gaseous TiO and VO would greatly heat the upper 
atmospheres of those planets already hot enough for these molecules to be present as vapor. 
Thermal inversions on transiting hot Jupiters were first seen by Richardson et al. (2007) 
for HD 209458b and Harrington et al. (2007) for HD 149026b. The observed flux ratio at 
8 /im for HD 149026b agreed only with models that included a thermal inversion (Fortney 
et al. 2006). Temperature inversions have since been confirmed by Spitzer observations of 
HD 209458b (Knutson et al. 2008a), XO-lb (Machalek et al. 2008), and TrES-4 (Knutson 
et al. 2009), all of which show distinctive flux ratios in IRAC bands that suggest inversions 
(Fortney et al. 2006; Burrows et al. 2007). More circumstantial evidence exists for HD 
179949b (Barnes et al. 2008). 

On the other hand TrES-1, the least irradiated planet with published Spitzer 
observations, does not appear to have a pronounced inversion (Burrows et al. 2008). Nor, 
seemingly, does HD 189733b, which is also modestly irradated (Charbonneau et al. 2008; 
Barman et al. 2008). One suggestion is that temperature inversions are triggered by 
irradiation reaching a critical level that is hot enough to evaporate TiO and VO from 
grains, as discussed by Burrows et al. (2007), Fortney et al. (2008), and Burrows et al. 
(2008). However, irradiation of XO-lb and HD 189733b is within uncertainties the same 
(Torres et al. 2008), which poses a challenge to the irradiation trigger. 

In the Solar System, stratospheric temperature inversions are often caused by 
absorption of UV light by gases or aerosols produced by photochemistry. Here we ask if 
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atmospheric chemistry might play a similar role in hot Jupiters. Speculation has tended to 
focus on sulfur-containing species (Tinetti 2008), as the reservoir species H 2 S is expected 
to be abundant (Visscher et al 2006) in these atmospheres and many of its breakdown 
products (S2, in particular) absorb violet and ultraviolet light. 

2. The Photochemical Model 

Previous photochemical modeling of hot Jupiters addressed the abundance of 
photochemical H (Liang et al 2003) and the absence of photochemical smogs (Liang et al 
2004). Liang et al (2003) focused on the high H/H 2 ratio that arises from H 2 photolysis. 
In their second paper, Liang et al (2004) argued that simple hydrocarbons would not 
condense to form photochemical smogs in hot solar composition atmospheres. Neither study 
considered sulfur. 

We have developed a new general purpose ID photochemical kinetics code applicable 
to hot extrasolar planets. The code is based on the sulfur photochemistry model for early 
Earth originally described by Kasting et al (1989) and Kasting (1990), and subsequently 
adapted by Zahnle et al (2006) and Claire et al (2006) to address sulfur photochemistry of 
Earth's atmosphere during the Archean, and by Zahnle et al. (2008) to address martian 
atmospheric chemistry. Steady state solutions are found by integrating the system through 
time using a fully implicit backward-difference method. 

Our chemical network has been upgraded from that used by Zahnle et al (1995) to 
address the chemistry generated when the fragments of Comet Shoemaker Levy 9 struck 
Jupiter. We have assembled a reasonably complete list of the reactions that can take place 
between the small molecules and free radicals that can be made from H, C, O, N, and 
S. The code solves 507 chemical reactions for 49 chemical species: H, H 2 0, OH, O, 2 , 
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CO, CO2, HCO, H2CO, C, CH, CH2, CH3, CH4, CH3O, C2, C2H, C2H2, C2H3, C2H4, C 
2 H 5 , C 2 H 6 , C 4 H, C4H2, CN, HCN, N, N 2 , NO, NH, NH 2 , NH 3 , NS, H 2 S, HS, S, S 2 , S 3 , 
S4, S 8 , SO, HSO, S0 2 , OCS, CS, HCS, H 2 CS, CS 2 , and H 2 . Reaction rates, when known, 
are selected from the publicly available NIST database (http://kinetics.nist.gov/kinetics). 
In order of decreasing priority, we choose between reported reaction rates according to 
relevant temperature range, newest review, newest experiment, and newest theory. Reverse 
reaction rates k r = K eq kf of two-body reactions are determined from the forward reaction 
rate kf and the equilibrium K eq = exp {(—AH + TAS) / RT} by using H°(T) and S°(T) 
as available (R is the gas constant). Rates are not available for all reactions, especially for 
reactions involving elemental sulfur. We will present a full listing of the chemical reactions 
important to sulfur in a more general followup study. 

Here we use simple descriptions of atmospheric properties. The background atmosphere 
is 84% H 2 and 16% He. We include Rayleigh scattering by H 2 (Dalgarno and Williams 
1962). For our base case we assume an isothermal atmosphere with T=1400 K; constant 
vertical eddy diffusivity K zz = 1 x 10 7 cm 2 /s; a surface gravity of 20 m/s 2 ; and insolation 
levels / by a solar twin that are 1000 x greater than at Earth. Metallicity proved to be 
the most interesting parameter and was varied —0.3 < [M/H] < 1.7. In these units, solar 
metallicity is [M/H] = 0, Jupiter's is [M/H] = 0.5, and Saturn's is [M/H] = 0.8. Short 
chemical lifetimes of S-containing species make our results insensitive to K zz . Model 
parameters are listed in Table 1. 

At the upper boundary we set a zero flux lid at 1 /xbar, with neither escape nor 
exogenous supply. For the lower boundary we use fixed equilibrium mixing ratios of the 
most abundant species at 1 bar of H 2 and temperature T (Lodders and Fegley 2002, 
Visscher et al 2006). For other species we force the mixing ratio at 1 bar to approach zero. 
We scale the lower boundary conditions such that the total mixing ratios of C, O, N, and S 
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all scale linearly with metallicity. 

Absorption by S 2 between 240 nm and ~360 nm from the ground state is analogous to 
the Schumann- Runge system in 2 (Okabe 1978). Strong, distinctive S 2 emission near 300 
nm was observed on Jupiter after the impact with Shoemaker-Levy 9 with Jupiter in 1994 
(Noll et al 1995). Subsequent thermochemical modeling showed that S 2 readily forms as a 
major product in a shock-heated (T > 1000 K) gas of either cometary or jovian composition 
(Zahnle et al 1995, Zahnle 1996). S 2 has also been seen in gases vented by volcanoes on Io 
(Spencer et al, 2000; Moses et al 2002). For S 2 , we use absorption cross sections at 1500 K 
computed by van der Heijden and van der Mullen (2001). 

The HS (mercapto) radical absorbs from its ground state at 324 nm (Okabe 1978). 
Visscher et al (2006) predicted that HS would be very abundant in equilibrium at hot 
Jupiter conditions. We find the same. We therefore calculated absorption cross sections 
of HS at four temperatures at 30 mbar pressure using literature values of the molecular 
properties. The ground X 2 Ii state has been well studied (Ram et al 1995) but the upper 
level A 2 S + is subject to strong pre-dissociation (Resende and Ornellas 2001, Wheeler et al 
1997, Schneider et al 1990, Henneker and Popkie 1971), and only the value of the rotational 
constant B and the spacing of the lowest vibrational energy levels have been well measured. 
Using these constants and a value for the electronic band oscillator strength of the 0-0 
transition derived from a study of HS in the solar spectra by Berdyugina and Livingston 
(2002), a line list was computed using the RLS code developed by R.N. Zare and D. 
Albritton (Zare et al 1973). This RLS code uses the molecular constants and band strengths 
to predict line positions and strengths by fitting to an RKR potential (Zare et al 1973). 
Other needed data — Franck-Condon factors, partition functions, etc. — were derived 
either from the cited literature, the program itself, or from Sauval and Tatum (1984) or 
Larsson (1983). The calculations were carried out for values of t>"(0-4) and v '(012). Because 
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the excited vibrational levels of the A 2 H + state are unstable with respect to predissociation, 
the corresponding optical transitions are likely to be broad and shallow, or even continuous. 
These uncertainties principally affect the absorption spectrum at wavelengths shorter than 
324 nm, which is in the range that is absorbed strongly by S2. Results are shown in Figure 
[TJ In the photochemical model we used only the 1500 K absorption coefficients. 

Other sulfur allotropes are better absorbers than S2 but less abundant. S3 absorbs 
strongly between 350 and 500 nm, and S4 absorbs between 450 nm and 600 nm, but more 
weakly (Billmers and Smith 1991). Unfortunately, the chemistries of S3 and S4 are very 
uncertain, and we have had to estimate the important reaction rates. In an earlier version 
of this study, we focused on the heats of formation, and we tentatively concluded that 
S3 heating would be important for metallicities [S/H] > 0.7. We have since learned that 
reactions of the form H + S n — > HS + S n _i, where n > 2, are strongly favored by entropy. 
The revised model predicts less S2 and much less S3, which reduces the importance of S3 
heating considerably. 

Sulfanes (H 2 S n , hydropolysulfides) will be present in cooler hot Jupiters. At low 
temperatures sulfanes absorb VUV between 260 nm and 330 nm (Steudel and Eckert 2003). 
Absorption may extend beyond 400 nm at higher temperatures as the ground state becomes 
vibrationally excited, as in HS, but to first approximation these wavelengths are covered by 
the more abundant S2 and HS. We have not included sulfanes in this study. 

3. Results 

Figure [2] shows how CO2 and the abundant S-containing species vary as a function of 
altitude. This particular case shows a hot Jupiter at 1400 K with a "planetary" metallicity 
of [M/H] = 0.7. Figure [2] is broadly representative of all our models with 1200 <T< 2000 
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K and —0.3 < [M/H] < 1.7. In particular, S 2 and HS show well-defined peaks at ~2 mbars 
that coincide with the altitude where H 2 S photolysis becomes important. At lower altitudes 
H 2 S is the main S-containing species, and at higher altitudes S is. It is also notable that the 
atmosphere becomes more oxidizing at higher altitudes where H 2 photolysis is important. 

Table 1 lists some key results pertinent to sulfur for several variations of basic model 
parameters. The models assume that K zz = 10 7 cm 2 /s and g = 2000 cm/s 2 unless otherwise 
noted. In this temperature range the models are insensitive to K zz (results not shown). 
Model G shows that, as expected, column densities vary inversely with g. 

Column densities of S 2 and HS are sensitive to metallicity. To first approximation, 
species with one metal atom, such as H 2 and H 2 S, increase linearly with metallicity, and 
species with two metal atoms, such as SO and S 2 , increase as the square of metallicity 
(VIsscher et al 2006). A slight complication is that CO and N 2 increase linearly with 
metallicity because these are the major reservoir species for C and N, respectively; hence 
C0 2 increases as the square of metallicity (as CO x O), rather than as the cube. 

The models are not sensitive to temperature and insolation over the parameter ranges 
(1200 < T < 2000 K and 1 < / < 1000) presented here. Insensitivity of the chemistry to 
T and / surprised us, and suggests that thermochemical equilibrium is more important for 
sulfur than photochemistry or kinetics. Minor differences are that HS is favored by higher 
temperatures and SO and S 2 are favored by high /. Not shown here is that the chemistry 
changes markedly for T < 1100 K: hydrocarbons, CS, and CS 2 become abundant, and the 
results become sensitive to K zz . Cooler atmospheres introduce a variety of new topics best 
left for another study. 

Carbon dioxide, a robust molecule and a potential observable, has been reported 
in HD 189733b by Swain et al (2009). C0 2 is generated from CO by reaction with OH 
radicals. The chief source of OH is the reaction of H 2 with atomic hydrogen; at high 



- 9- 



altitudes UV photolysis of H2O is also important. We find that CO2 mixing ratios range 
from 1.6 x 10~ 8 to 1.6 x 10~ 4 for —0.3 < [M/H] < 1.7, scaling as the square of metallicity. 
Table 1 lists computed CO2 mixing ratios in the models discussed here. These results 
are insensitive to insolation, vertical mixing, temperature between 1200 K and 2000 K, 
and gravity. The CO2/CO ratio is nearly independent of pressure, as seen in Figure [21 
Pressure independence is expected because the controlling reactions, C0 2 +H «-> CO + 
OH and H+H 2 <-> H 2 + OH, and the controlling equilibrium, C0 2 +H 2 <-> CO + H 2 0, 
all leave the total pressure unchanged. (At very high altitudes photochemistry alters the 
CO2/CO ratio.) The computed CO2 abundances are in good agreement with the reported 
observation of CO2 at the ppmv level in HD 189733b (Swain et al. 2009). The sensitivity 
of CO2 to metallicity and insensitivity to other atmospheric parameters makes CO2 a good 
probe of planetary metallicity, as pointed out by Lodders and Fegley (2002). 

3.1. Optical depth and stratospheric heating 

Figure [3] shows the pressure levels where the solar and planetary metallicity atmospheres 
of Models A, M, and MM become optically thick. Opacity is dominated by HS, with some 
contribution by S2 at wavelengths shorter than 300 nm. The twin peaks between 300 nm 
and 320 nm may be fictitious, but the peak at 324 nm could prove diagnostic of HS. A solar 
and a K0V stellar spectrum are shown for comparison. 

Figure H] shows the magnitude of stratospheric heating and the pressure level where 
the heating occurs for a solar-twin primary at 0.032 AU (/ = 1000) for 3 metallicities 
(Models A, M, and MM). Radiative heating is dominated by HS, and is nearly saturated 
through the stratosphere for all these models (see also Table 1). By contrast, peak heating 
at ~ 100 /ibars takes place where SO and SO2 are significant. The sensitivity of SO and 
SO2 to metallicity is reflected in greater heating rates at ~100 /ibars. 
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Cumulative stratospheric heating rates for these models are listed in Table 1. For a 
solar-twin at 0.032 AU, cumulative heating above 1 mbar is typically 4 x 10 4 W/m 2 and 
above 0.1 bars is typically 8 x 10 4 W/m 2 , i.e., about half the energy is absorbed in the lower 
stratosphere. Burrows et al. (2008) modeled hot stratospheres by adding an unknown gray 
absorber. They found that gray cross-sections of 0.05 — 0.6 cm 2 /g, averaged over 430 to 
1000 nm for altitudes above 0.03 bars, could produce the observed heating. Heating profiles 
using gray opacities in this range are plotted for comparison on Fig H] for the same planet 
and star. The gray opacities produce more heating in total (indeed, the stratospheres in 
both these models are optically thick), and more heating at low altitudes, but at higher 
altitudes sulfur generates heating at levels quite similar to what Burrows et al find useful. 

4. Conclusions 

We develop a new ID photochemical model for stratospheric modeling of hydrogen-rich 
atmospheres of warm or hot exoplanets. This model is applicable to any H-rich planet 
subject to high insolation, including hot Neptunes, superearths, and waterworlds. Here we 
apply the model to sulfur chemistry, stratospheric heating, and CO2 abundance. 

We find that hot stratospheres of hot Jupiters could be explained by absorption of 
UV and violet visible light by HS and S2, two highly reactive species that are generated 
chemically from H 2 S. For a hot Jupiter orbiting a solar-twin at 0.032 AU, for a wide range 
of possible planetary compositions, HS and S 2 together absorb 4 x 10 4 W/m 2 at altitudes 
above 1 mbar and another 4 x 10 4 W/m 2 at altitudes between 1 mbar and 0.1 bar. This 
level of heating approaches what Fortney et al (2006) and Burrows et al (2008) use in their 
most successful LTE spectral models. Non-LTE mechanisms may improve the agreement, 
because LTE models systematically overestimate radiative cooling and thus underestimate 
the temperature. Chemiluminescence by H 2 0, formed by the exothermic reaction of 
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OH+H 2 , might also be expected. 

Although our computed HS and S 2 column densities increase with metallicity, optically 
thick columns are predicted for all plausible atmospheric compositions, which means that 
millibar-level temperature inversions are expected to be commonplace. The distinctive 
interaction of S 2 and HS with near ultraviolet light could make these species detectable in 
transit by the refurbished HST; there is evidence for a blue absorber in legacy HST data of 
HD 209458b (Sing et al 2008). 

On the other hand, sulfur does not give an easy answer to why some hot Jupiters 
have superheated stratospheres, and others not. In an earlier draft of this study, we 
speculated that S3 — which is very sensitive to metallicity — might be part of the explanation. 
This no longer appears likely. We have since developed a better understanding of HS's 
opacity, which turns out to be considerable. We no longer see a strong connection between 
metallicity and radiative heating, save at very low pressures (< 100/zbars) where SO and 
SO2 become important. It now seems that sulfur chemistry by itself is unlikely to explain 
differences between planets, although planetary metallicity may still be key. 

Heating by sulfur compounds does not preclude heating by TiO and VO on hotter 
planets. Sulfur species provide considerable heating from below 1000 K to above 2000 K, 
but they do not provide the spectral coverage at visible wavelengths that TiO and VO 
provide. For TiO and VO to be abundant enough to explain stratospheric heating, the 
temperature needs to be very high, in excess of 2000 K, and not just in the stratosphere 
but also at deeper levels in the planet where these two refractory oxides would otherwise 
be cold-trapped in silicate clouds. OGLE-TR-56b (Sing and Lopez-Morales 2009) seems to 
meet the TiO-VO threshold. 

C0 2 is generated by the reaction of CO with OH and destroyed by the reverse 
(endothermic) reaction with H, CO + OH <-> CO2 + H. At low altitudes OH is generated 



- 12 - 



by the reaction of H 2 with atomic H, supplemented at high altitudes by UV photolysis 
of H 2 0. As both the major source and major sink of CO2 are proportional to atomic 
hydrogen densities, the kinetic inhibition against hydrogen recombination does not disturb 
CCVs thermochemical equilibrium. We find that CO2 mixing ratios vary quadratically with 
metallicity from 1.6 x 1(T 8 to 1.6 x 1(T 4 for < [M/H] < 0.7. This result is insensitive to 
insolation, vertical mixing, temperature (for 1200 < T < 2000 K), and gravity. Because 
the reactions that form and destroy C0 2 leave the total number of molecules unchanged, 
the C0 2 /CO ratio is also pressure independent. The computed C0 2 abundances are in 
good agreement with the observation of C0 2 at the ppmv level in HD 189733b (Swain 
et al 2009). Therefore we confirm Lodders and Fegley's (2002) suggestion that C0 2 is a 
promising probe of planetary metallicity. 
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10 21 


26 


1.7 


1.6 


8.4 x 


10 4 


4.1 x 10 4 


MM 


1.4 


1000 


1400 


6.0 x 


10 21 


5.4 x 


10 21 


100 


27 


41 


1.1 X 


10 5 


5.1 x 10 4 


H 


0.7 


1000 


1600 


2.2 x 


10 20 


2.3 x 


10 21 


43 


1.2 


1.4 


9.0 x 


10 4 


4.3 x 10 4 


HH 


0.7 


1000 


1800 


1.9 x 


10 20 


4.0 x 


10 21 


52 


1.3 


1.5 


9.5 x 


10 4 


4.4 x 10 4 


HHH 


0.7 


1000 


2000 


1.3 x 


10 20 


5.0 x 


10 21 


41 


1.4 


1.3 


9.6 x 


10 4 


4.1 x 10 4 


C 


0.7 


1000 


1200 


1.6 x 


10 20 


2.5 x 


10 20 


11 


1.9 


1.9 


7.5 x 


10 4 


3.7 x 10 4 


G 


0.7 


1000 


1400 


4.3 x 


10 20 


2.0 x 


10 21 


32 


1.3 


1.6 


9.3 x 


10 4 


4.7 x 10 4 


I 


0.7 


200 


1400 


2.1 x 


10 20 


1.0 x 


10 21 


37 


0.9 


1.6 


1.7 x 


10 4 


8.8 x 10 3 


ssc 


0.7 


1 


1200 


1.1 X 


10 20 


2.4 x 


10 20 


16 


0.06 


1.9 




72 


44 



a - Metallicity. This notation means that the planet is 10[ M / H 1 richer in C, S, N, and O than the Sun. 
b - Insolation. I = 1000 corresponds to a solar twin primary at 0.032 AU. 
c - Column densities above 1 bar. 
d - Mixing ratio in ppmv at 1 mbar. 

e - Total atmospheric heating [W/m 2 ] above 0.1 bar for a solar twin source. 
/ - Total atmospheric heating [W/m 2 ] above 1 mbar for a solar twin source. 
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Fig. 1. — Theoretical absorption cross sections of HS radicals at near UV, violet and indigo 
wavelengths at four temperatures at 30 mbar pressure. Cross sections were computed from 
the lowest five vibrational levels of the ground electronic state X 2 U to the lowest three 
vibrational levels of the upper level A 2 Y> + . The excited vibrational levels of A 2 S + are 
strongly predissociating, which suggests that absorption at wavelengths shorter than 324 
nm is probably continuous rather than allocated into the well-defined bands shown here. 
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Volume Mixing Ratio 

Fig. 2. — Important sulfur species, CO, and CO2 in the atmosphere of a hot Jupiter with a 
"planetary" metallicity of [M/H] = 0.7. The atmosphere is assumed isothermal at 1400 K and 
insolated 1000 x more strongly than Earth. Other model M parameters are listed in Table 1. The 
prominent transition at ~2 mbar — where the S2 mixing ratio peaks — is associated with photolysis 
of H2S. The bump in CO2 at 6 /ibars is attributable to photochemistry. Abundance profiles in the 
1400 K atmosphere are generally representative of atmospheres with 1200 <T< 2000 K. 
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Fig. 3. — Pressure levels of the r = 1 surface as a function of wavelength for three metal- 
licities, [S/H] = 0, 0.7, and 1.4. These metallicities correspond to models A, M, and MM 
of Table [5j Absorption between 250 and 300 nm is mostly by S2 and absorption between 
300 and 460 nm is by HS. Structure blueward of 324 nm is associated with transitions to 
predisso dating states and is probably fictitious. The r = 1 surface of a pure H 2 Rayleigh 
scattering atmosphere and two incident stellar spectra, one for the Sun and another for a 
generic K0V dwarf, are shown for comparison. 
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Fig. 4. — Radiative heating at different altitudes for three metallicities, [S/H] = 0, 0.7, and 
1.4. These correspond to models A, M, and MM of Table Heating rates are given in 
W/kg, which emphasizes the potential impact on temperature. Heating peaks at 100 //bars 
but extends through the stratosphere. Heating with constant gray opacities of 0.05 and 0.6 
cm 2 /g for 430 < A < 1000 nm is shown for comparison. 



